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ABSTRACT 

We present near and mid-IR observations of a sample of Seyfert II galaxies drawn from the 12/xm 
Galaxy Sample. The sample was observed in the J, H, K, L, M and N bands. Galaxy Surface 
Brightness Profiles are modeled using nuclear, bulge, bar (when necessary) and disk components. 
To check the reliability of our findings the procedure was tested using Spitzer observations of M31. 
Nuclear Spectral Energy Distributions (SEDs) are determined for 34 objects, and optical spectra 
are presented for 38, including analysis of their stellar populations using the STARLIGHT spectral 
synthesis code. Emission line diagnostic-diagrams are used to discriminate between genuine AGN and 
HII nuclei. Combining our observations with those found in the literature, we have a total of 40 SEDs. 
It is found that about 40% of the SEDs are characterized by an upturn in the near-IR, which we have 
quantified as a NIR slope a < 1 for an SED characterized as A/a c>c A". Three objects with an HII 
nucleus and two Seyfert nuclei with strong contamination from a circumnuclear starburst, also show 
an upturn. For genuine AGN this component could be explained as emission from the accretion disk, 
a jet, or from a very hot dust component leaking from the central region through a clumpy obscuring 
structure. The presence of a very compact nuclear starburst as the origin for this NIR excess emission 
is not favored by our spectroscopic data for these objects. 
Subject headings: 
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1. INTRODUCTION 
1.1. Unified Model 

The presence of well-collimated radio jets in radio loud 
Active Galactic Nuclei (AGN), and the expectation of 
an axis of symmetry defined by the geometrically thin 
accretion disk, argue for strong anisotropics in the cen- 
tral engines of AGN. This implies that the appearance 
of an AGN will depend on the observer's position rela- 
tive to the system's axis. Following this reasoning, and 
to explain a wealth of observational evidence, the Uni- 
fied Model postulates that Type I and Type II AGN are 
the same kind of object, but that they appear distinct 
because of the difference in angle between the central 
engine's axis and our line of sight. 

Support for this model has been particularly strong in 
the case of Seyfert galaxies. Optical spectra of Seyfert I 
objects show Broad Emission Line Regions (BLR); those 
of Seyfert II objects do not. Explaining this as an orien- 
tation effect, wit h obscuration of th e central source, goes 
back to at least lOsterbrockl (|1978D . who suggested that 
both types are physically the same, but that in Seyfert 
II galaxies the BLR is blocked from our view, probably 
by dust in a toroidal structure surrounding the central 



source. Below we summarize some of the work that fol- 
lowed and supported this thesis. 

Evidence for obscuration of Type II objects comes from 
the excess absorption in X-ray, UV and optical wave- 
lengths , when observation s are compared with Type I 
nuclei (iLawrence fc Elvii 119821: iMass-Hesse et all (l993l: 



Turner et all 119971: IMalkan et al.l 119981: IRisaliti et alj 
1999D. Receiitly,[Sh i et al. I (120060 have shown a clear cor- 
relation between the column densities derived from X-ray 
data and the strength of the 9.7/im silicate feature: Type 
II systems with large hydrogen columns show large ab- 
sorption features, while objects classified as Type I often 
show the feature in emission. 

In those AGN that are radio-loud sources, the radio 
structure appears to be aligned with the torus axis (Wil- 
son & Tsvetanov, 1994; Nagar et al. 1999; Barbosa 
et al. 2009). This anisotropy is also observed in the 
cone-like structures of the Narrow Line Region, where 
the tor us coUimates io nizing radiation fr om the central 
source (|Poggelll988allbl : IFalcke et al.l[T99l . 

In the IR, detection of broad components from hy- 
drogen recombination lines such as Pa/3 and Br7 in ob- 
jects optically classified as Seyfert II galaxies is consis- 
tent with extinction along a line-of-sight that grazes the 



2 



Videla et al. 



torus surface, sufficient to suppress the BLR emission 
at optical wavelengths, but optically thin in the infrared 
(|Ruiz et al.llT99l IVeilleux et al.lfT997l[l999l) . 

More direct support for unification conies from the 
HST, which has imaged an extended nuc lear structure of 
~90 pc in NGC4261 (Jaf fe et aLl l TOO g). while interfer- 
ometry observations are starting to shed light on the in- 
ner obscuring torus (Wittkowski et al. 1998, 2004; Wein- 
berger et al. 1999; Swain et al. 2003; Jaffe et al. 2004; 
Weigeh et al. 2004; Meisenheimer et al. 2007; Tristram et 
al. 2007, 2009; Beckert et al. 2008; Kishimoto et al. 2009, 
2011; Raban et al. 2009). 

There is also evidence for structure in the ob- 

scuring material. Changes in broad line profiles 

(jPenston fc Perej 119841) and inferred column densities 
imply the presence of clumpy structures, leading to vari- 
able obscuration of the central engine as th ey move across 
the line of sight (.Risaliti c t al. 2011. 2009^ iPuccetti etall 
I2007t lElvis et al.l I2004D falthough the column density 
variations are more likely to correspond to obscuration 
by BLR clouds than the torus). 

Finally, the most famous, and arguably most convinc- 
ing evidence in support of the Unified Model comes 
from optical spectropo larimetry of Type II objects. 
lAntonucci fc Milleil (|1985i ) found scattered emission from 
the BLR in polarized flux from NGC 1068, strongly 
supporting the idea of unification by orientation, and 
so implying the existence of an obscuring structure . 
Their work was e xtended by (jYoung et al.lll995l 119961 : 
iMoran et al1[2000l ). 

1.2. Spectral Energy Distribution 

The Spectral Energy Distribution (SED) is an impor- 
tant tool for studying the central engine in AGN. In most 
luminous and intermediate luminosity AGN it can be 
characterized by two bumps — one peaking in the UV 
and the other in the mid-IR — and a n infiection point 
between, at ~ 1.5 /Ltm (|Elvis et al.lll994[ ). The UV feature 
(the "Big Blue Bump" ) is normally identified with emis- 
sion from the central accretion disk, while the IR flux 
is assumed to originate in the dusty obscuring structure, 
the torus. The latter has a broad SED, with a peak some- 
where between 5 — 30/im, decreasing slowly towards the 
far-IR. The emission mechanism is believed to be thermal 
reprocessing of central, optical-UV emission by dust at a 
wide range of temperatures — the highest 1000—1500 
K) corresponding to dust closest to the central engine. 

Obtaining the SED of the torus in AGN is difficult for 
a number of reasons. 

First, the SED inflection point matches the peak in 
starlight emission from the host galaxy. In intermediate 
luminosity Seyfert galaxies the host galaxy contribution 
is often comparable to the emission from the nucleus, 
flattening the optical-IR region of the SED and lowering 
the relative strength of the UV and IR bumps. To study 
the nuclear SED, galactic starlight must be removed. 

Second, until recently, only small size IR arrays have 
been available, making it impossible to image the com- 
plete host galaxy for nearby sources, compromising the 
modeling of the surface brightness distribution. 

Third, observations at A > lOyitm can only be achieved 
from space, which are characterized by lower angular res- 
olution. 

Fourth, and finally, observations obtained at different 



epochs can be affected by variability. Fortunately, the 
torus region is large, damping most variations of the cen- 
tral source, so this is not expected to be significant. 

The number of compiled AGN SEDs in the literature 
increases constantly, but many of them incorporate data 
gathered by different groups using different instruments; 
this often implies data of different resolutions, analyzed 
in different ways. To properly study the emission from 
the dusty torus in AGN, and to accurately infer physical 
and geometrical characteristics, this paper describes a 
more homogenous dataset. We have constructed nuclear 
IR SEDs of a sample of 48 Type II Seyfert galaxies using 
near- and mid~IR high-resolution ground-based images 
(described here), and mode l ed them in clu mpy using the 
approach of iNenkova et al.l (l2008allbl. I2002D (described in 
a companion paper. (jLira et al.l 120121 )). The main con- 
tribution of this work is the number of objects we study, 
and the detailed treatment of the nuclear SED construc- 
tion. 

In [J2] we present the sample, observations and data 
reduction, including the optical spectroscopy and the 
stellar population analysis. fj3] describes, in detail, the 
path from IR images to the final nuclear IR SEDs of our 
sample, including the analysis of M31 to validate our 
methodology. Section ^ discusses the results and possi- 
ble interpretations. A summary is presented in iJ5] 

2. OBSERVATIONS 

The sample used in this work contains all Seyfert 
II galaxies from the Extended 12/im Galaxy Sample 
(jRush et al.lll993l ) that are located in the southern hemi- 
sphere. The 48 objects are described in Table [H 

The 12/im Galaxy Sample has several advantages over 
other catalogs. Most importantly, it is selected in 
the mid-infrared (MIR), giving a more representative 
population of nearby (z<0.07) AGNs. In comparison, 
optically-selected samples can miss obscured systems 
and are more susceptible to the effects of host galaxy 
reddening (e.g., Huchra & Burg, 1992; Maiolino & Rieke, 
1995). This is important because a wide range of nuclear 
obscuration (hydrogen columns of ^ 10^^ — 10^^ cm~^) 
can provide more general tests of nuclear emission mod- 
els. In addition, the 12/im Sample includes elliptical, 
lenticular and spiral galaxies (which helps avoid system- 
atic errors in the decomposition process of the surface 
brightness profiles) and a wide range of galaxy inclina- 
tions (which allows the importance of galactic, as well as 
nuclear, absorption to be assessed). 

However, the sample may be biased towards Seyfert II 
sources with strong star formation (Buchanan et al. 2006) 
— the stellar IR emission can compensate for a weak 
AGN component, increasing the likelihood of inclusion in 
a flux limited sample. This might explain the prevalence 
of star formation in Seyfert II galaxies when compared 
with Seyfert I nuclei (e.g., Maiolino et al. 1995). 

Our analysis of the central region of galaxies is built 
on extensive studies of the 12/im Galaxy Sample at many 
wavelengths, from radio to X-rays (Spinoglio et al. 1995, 
2002; Rush et al. 1996b; Hunt et al. 1999ab; Bassani et 
al. 1999; Thean et al. 2000, 2001; Imanishi, 2003, 2004; 
Gorjian et al. 2004; Strong et al. 2004; Buchanan et 
al. 2006; Tommasin et al. 2008, 2010; Wu et al. 2009; 
Gallimore et al. 2010; Baum et al. 2010), including spec- 
tropolarimetric observations (Tran 2001). 
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TABLE 1 

Observed sample of 48 Seyfert 2 galaxies, indicating the detector used in each infrared observation. 





Hubble Type 


Redshift 






RA 




DEC 


JHK bands 


LM bands 


N band 






(■31 


(4.) 




(5) 




(6) 




\°) 




NGC34, MRK938 


iSc 
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02 
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- 
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99 
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02 
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40 


27 
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- 
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- 
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02 


00 
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iU 


85 
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0.64 
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UU 


oU 


69 
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SAO-I- 
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03 


HQ 
UO 




11 
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- 
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03 
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11 
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90 
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ISA-2 


- 


- 
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03 


OA 
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72 
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ISA-2 


- 
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NGC 1386 
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03 


OD 
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40 
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ISA-2 


- 


- 


F 03362-1642 
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03 


00 


04 


54 
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SO 
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04 
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04 


96 
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ISA-2 


- 
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SAB(r)c 
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04 


^0 


/ 


14 
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ISA-2 


- 
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ES0 33-G2 
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04 


00 


00 
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21 


12 
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NGC 4388 
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25 


46 


75 
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12 


31 


59 


16 
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0.12 


12 


40 


52 


88 
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13 


06 
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Note. — In column (2), i indicates Compton Thick (CT) sources. In column (4) we give the galaxy inclination obtained from the mean 
ellipticity (e = 1 — b/a) of the last fitted isophotcs in the near-IR bands ( ^3.11 1: 2 nuc indicates the detection of 2 nuclei in the center of 
the galaxy, close enough to prevent reliable SBP fitting of 2 separated galaxies; lost indicates that the observations of the galaxy were lost, 
either beca use of t he size of the galaxy compared to the detector (bad sky subtraction) , or because of non-photometric observing conditions 
(NGC 424, ^2.2.11 1. In column (7), V indicates observations made in visitor mode and S indicates observatio ns made in service m o de. I n 
column (9), G04 indicates observations made by Gorjian et al. (2004); AH03 indicates SEDs determined by lAlonso-Herrero et ahl II2003I) ; 
NED indicates values obtained from the NASA/IPAC Extragalactic Database" 

"The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Administration 

The initial classification of objects in the sample, 
by Rush et al. (1993), was based on existing cat- 
alogs of active galaxies (jVeron-Cettv fc VeronI I199TI : 
IHewitt fc Burbidgel[T99l . Since then, as better opti- 
cal spectra became available, some objects have been re- 
classified. To verify the classifications, and to study stel- 
lar populations, we obtained high quality optical spec- 
troscopy for 38 objects. The analysis of these data is 
presented next. 



2.1. Optical Spectroscopy and Analysis 

Long-slit optical spectroscopy was obtained with the 
R-C Spectrograph on the Blanco 4m telescope at CTIO, 
in August 2007 and February 2008. We used the KPGL3- 
1 grating and a slit width of 1". The slit was positioned 
at the parallactic angle to minimize light losses. In both 
runs the seeing varied between ~0.8 and ^ 1.2 arc-sec. 
The February 2008 run was photometric, while clouds 
were presented during the August 2007 run. However, 
absolute flux calibration is not necessary to characterize 
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Fig. 1. — Diagnostic diagrams for emission line spectra. The 
theoretical divisions between the star-forming, AGN and LINER 
classes, as parameterized by Kcwlcv et al. (2001), are shown with 
thick gray lines in panels (a), (b), and (c). Dotted lines indicate 
±0.1 dex. Nuclei that fall in the HII region are labeled. Also, AGN 
nuclei are shown with an extra circle or square if they have a NIR 
slope grater or smaller than 1, respectively (see §4). 



the nebular and stellar components of the nuclear spec- 
tra. 

The data were reduced using IRAF^ software. Stan- 
dard reduction procedures (bias subtraction, flat-fielding, 
slit illumination correction) were followed. Extrac- 
tion of the galactic spectra used a window of width 
~2.2". Wavelength calibration used comparison copper- 
neon and copper-argon lamps. The spectra were 
corrected for the f oreground extinction tabulated by 
ISchlegel et aT (199^ using the empirical extinction func- 
tion of Cardclh et al. (1989). 

We determined the stellar component of the spectra us- 
in g the STARLIGHT spectral synth esis code developed 
bv lCid Fernandes et al.l ([200l [2001 . A total of 45 dif- 
ferent stellar population templates were used, covering 
an age range of ^ 10^ — 10^° years and a metalicity range 
of 0.2 — 2.5 X Zq. The best fit spectrum was found by 
minimizing \ ^ after masking all em is sion li nes. 

Following iCid Fernandes et al.l (|200l [2005) we 
grouped output from the population synthesis into 3 
age groups: a Young component (Y: t < 10^ yr); an 
Intermediate component (I: 10^ < t < 10^ yr); and an 
Old component (O: t > 2.5 x 10^ yr). The results are 
presented in Table [2] 

The best-fit stellar models were subtracted from the 
observed data to give pure emission line spectra. This 
reveals weak emission lines present in targets with sig- 
nificant intermediate age stellar populations. Table [2] 
contains the measured emission line flux ratios. 

Diagnostic diagrams were constructed using emission 
line fluxes and are shown in Fig. [T] The diagrams al- 
low us to distinguish different ionizing mechanisms of 
the emission nebula using line ratios that are sensitive 
to changes in the spectra l shape of the ionizing radia- 
tion (iBaldwin et all [TOSlt iVeilleux fc Osterbrockl [l987l : 
IKewlev et al.l 120011 ). In this way we can test the identi- 
fication of the sources as Seyfert nuclei. 

We will discuss possible contaminants in Fig.[T]by type. 
First, we argue that there is little evidence for LINERs 
in our sample. 

In panel (a), there are no data points in the LINER 
region. The border between the HII and LINER regions 
is defined to include all star-forming models (it is cal- 
culated from the highest possible line ratios), but the 
[Nil] A6584/Ha ratio for star-f ormin g regions and LIN- 
ERs can overlap (jKewlev et al.ll2001| ) , so LINERs can be 
located to the left of the border. Also, theoretical models 
of shock-excitation give line ratios in the LINER regio n 
only when no precursor is observed ([Alien et al.ll2008D . 
For both these reasons panel (a) does not discriminate 
these two classes of object well. 

In panel (b), variations in density can introduce un- 
certainties in the [SII] AA6717,6730/Hq! ratio, so objects 
should not be classified as a LI NER based purely on this 
diagnostic ([Kewlev et al.ll200lD . 

Given the above, and the single LINER candidate in 
panel (c) which gives the most reliable test, none of our 
sources can be firmly classified as a LINER nucleus. 

Second, four sources systematically appear in the HII 



^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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TABLE 2 

Population synthesis analysis and measured emission line fluxes and flux ratios. 



Name 


Y{%) 
(2) 


I(%) 
(3) 


0(%) 
(4) 


/(4861) 
(5) 


/(5007) 
(6) 


(5007/4861) 
(7) 


(6583/6563) 
(8) 


(6300/6563) 
(9) 


(6716/6563) 
(10) 


(6563/4861) 
(11) 




46.8 


34.6 


18.6 


0.9 


1.6: 


0.2: 


-0.02 


-1.1 




13.7 


b uuiy»-7yzD 


0.0 


0.0 


100.0 


0.5 


3.9 












r UU5zl-7Uo4 


0.0 


10.6 


89.4 


0.8 


6.8 


0.9 


-0.05 


-0.9 


-0.6 


8.3: 


T71 c r\ c -1 1 T/^ 1 
h/bU o41-lGlz 


53.6 


0.0 


46.4 


0.7 


5.8 


0.9 


-0.2 


-1.2 


-0.6 


6.1: 


ivT /1 y1 


23.9 


0.0 


76.1 


10.0 


54.5 








~z 




1^ Ui4ro-U/4U 


53.6 


0.0 


46.4 


1.4 


7.5 


0.7 


-0.2 


-0.9 


-0.7 


7.1 
























■\Tr^ 1 nor/ 

INCjO 1U97 




















"~ 


JNCjO 1125 


50.3 


3.9 


45.8 


2.4 


14.7 


0.8 


-0.09 


-0.9 


-0.3 


4.8 


IVT 1 1 A A 

In LrO 1144 


0.0 


0.0 


100.0 


0.6 


6.4 


77 


77 








MULj -z-o-oy 


21.8 


0.0 


78.2 


1.8 


26.8 


1.2 


-0.09 


-0.8 


-U.4 


4.Z 


In LtU liy4 


26.7 


0.0 


73.3 


0.3 


2.6 


0.9 


-0.3 


-0.4 


-0.3 


D.D 






~ 




"~ 


"7 


7Z 




~ 








4.80 


15.9 


79.3 


2.6 


25.5 


1.0 


-0.2 


-1.0 


-U.4 


3.8 




"~ 






"~ 


~ 




77 


~ 




"~ 


r Uv5v>DZ-1d4Z 


51.6 


0.0 


48.4 


0.7 


4.5 


0.8 


-0.09 


-1.1 


~z 


5.3 


1 U4oo0-Uoz5 


45.8 


0.0 


54.2 


0.5 


1.8 


0.5 


-0.1 


-0.8 


-0.02 


A n. 

4.U: 


■\T 1 c 


25.2 


0.0 


74.8 


1.1 


7.9 


0.9 


0.2 


-0.7 


-0.1 


5.7: 




32.3 


1.2 


66.5 


1.3 


VIA 


1.0 


0.07 


-1.0 


-0.2 


5.0: 


r Uoloy-ZOz4 


1.9 


26.7 


71.4 


1.9 


11.0 














75.6 


0.0 


24.4 


2.7 


15.1 


0.7 


-0.3 


-1.0 


-0.3 


3.6: 


"N/Tr^/^ 1 n on oq 


28.2 


10.7 


61.1 


0.8 


0.9: 


0.02; 


-0.2 


-1.2 


~z 


8.1 


IVir^/^ Of? en 

INCjO oddU 


3.3 


0.0 


96.7 


1.3 


3.5 


0.4 


-0.2 


-1.1 


-0.6 


3.3: 


ivT ^ Q Q 


35.2 


0.0 


64.8 


4.9 


45.2 


1.0 


-0.2 


-0.8 


-0.1 


A n 
4.U 


In LtO 4oU1 


7.8 


0.0 


92.2 


0.9 


3.6: 


0.6; 


0.3 


-0.8 


0.04 


6.1. 


T'OT 10QQ QCIA 

I UL, 1Zv>o-od4 


26.2 


20.5 


53.3 


8.1 


49.8 


0.8 


-0.08 


-0.8 


n A 
-U.4 


A n 
4.U 


1S.T ^ n ^ -i 


U.U 


U.U 


1 nn n 
lUU.U 


Q 
Z.O 


oi.y 


1 n 
i.U 


1 -1 
-i.i 


-U. / 


n 1 
-U.l 


A a 
4.0 


i\ LtU 4yDo 


28.7 


6.2 


65.1 


2.4 


26.0 


1.0 


0.07 


-0.8 


-0.3 


A Q 
4.0 


iVlULx -o-o4-d4 


0.0 


22.0 


78.0 


16.3 














IN Kjyj OloO 


OvS.y 


1 ^1 n 
iD.y 


on 


4.0 


OQ Q 


0. / 


n no 


-1.^5 


-U.O 


0.0 


MHlv 4Doa 


60.1 


23.7 


16.2 


1.1 


2.5 


0.4 


-0.2 


-0.9 


-0.2 


3.2 


JVIHK 4Dob 


63.6 


19.1 


17.3 


9.2 


67.0 








"T 




INLtO ooUd 


53.1 


0.0 


46.9 


9.2 


59.7 


0.9 


0.0 


-0.8 


-0.1 


7~A 

4.4 


JNLtO 5950 


66.6 


12.2 


21.2 


8.0 


0.3: 


-0.7: 


-0.4 




-0.6 


4.2 


T\ T .in J 

MCG -2-40-4 


27.3 


4.2 


68.5 


1.6 


8.8 












1 Io4ou-Uo44 


11.2 


13.8 


75.0 


1.7 


65.0 


1.1 


-0.1 


-0.9 


-0.7 


4.0: 


r iyzo4- / z4o 


18.9 


0.0 


81.1 


4.7 


4.3 












NGC 6810 


44.2 


2.1 


53.7 


0.6 


1.1: 


-0.3; 


-0.2 


-1.3 


-0.5 


7.4: 


NGC 6890 


19.2 


6.2 


74.5 


2.5 


15.0 


1.2 


0.09 


-0.9 


-0.5 


5.0 


IC 5063 








1.0 














MRK897 


34.3 


38.2 


27.5 




1.5 


-0.3 


-0.3 


-1.5 


-0.6 


4.3 


NGC 7130 


65.8 


12.1 


22.1 


3.0 


27.1 


0.8 


0.08 


-1.1 


-0.6 


7.9 


NGC 7172 








4.8 














MCG -3-58-7 


33.6 


27.8 


38.7 




18.6 


1.0 


0.09 


-1.0 


-0.4 


4.7 


NGC 7582 








5.9 














NGC 7590 


23.0 


8.0 


69.0 




3.3 


0.8 


-0.04 


-1.0 


-0.2 


5.9 


NGC 7674 








0.5 














CGCG 381-051 























Note. — Columns (5) and (6) are in units of lO"^** ergs/s/cm^. Unreliable flux measurements (blending and calibration issues) are 
labeled with ';'. Columns (7) to (10) report the logarithm of the flux ratios. Column (11) is the Balmer decrement. *; Spectra for both 
nuclei in MRK 463 were obtained. 



section of the diagrams, consistent with the classification 
of Gallimore et al. (2010). These are: NGC 5953 (but see 
next), MCG-HO-29-23, NGC 6810, and MRK 897. Also, 
NGC 34 appears systematically close to the starburst re- 
gion. 

NG C 5953 was classified as an AGN bv lRafanelli et afl 
(|199Cir ) after making a careful identification of the nuclear 
source, which can be confused with a bright star-forming 
region located 3" to the West. Unfortunately we centered 
our slit on the star-formation knot and the slit orienta- 
tion did not include the nucleus, so our optical spectrum 
is of the wrong object. The IR images, however, are dom- 
inated by the nuclear region — the star-forming region 
can also be seen, but is much less prominent. The de- 
termined SED, however, suggests that contamination by 



this circumnuclear starburst or a compact nuclear star- 
burst is present (§4.1). 

One of the most reliable signatures of an active nucleus 
is the detection of luminous, hard X-ray emission. We 
checked the literature for X-ray data on MCG +0-29-23, 
NGC 6810, and MRK 89 7, but found littl e support for 
classification as Seyferts. iStricklaiidI ()2007D presents the 
discovery of a starburst driven wind in NGC 6810 from 
XMM-Newton observations, and explains that the claim 
of a hard ^ 10^^ ergs/s source in this galaxy was 
due to the misidcntification of a HEAOl source. Upper 
limits of Lx < 1 X 10"^^ and 6 x 10^^ erg s/s in the soft 
band were reported by iRush et all (|1996l ) for MRK 897 
and MCG +0-29-23. 

Given the above, NGC 5953 seems to harbor a genuine 
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TABLE 3 

Summary of infrared observing runs. 



Period 


Telescope 


Detector 


Mode 


Bands 


2002B 


ANTU (UTl) 


ISAAC 


Visitor 


JHK 


2002B 


ANTU (UTl) 


ISAAC 


Service 


LM 


2003A 


ANTU (UTl) 


ISAAC 


Service 


LM 


2003B 


NTT 


Son 


Service 


JHK 


2004A 


Gemini S 


T-ReCS 


Service 


N 


2004B 


NTT 


Son 


Visitor 


JHK 



TABLE 4 
Summary of the seeing conditions. 



Band 


Mean 


cr 


Median 


Mean 


cr 


Median 




["] 


["] 


["] 


[pc] 


[pc] 


[pc] 


J 


0.7 


0.1 


0.7 


325 


317 


224 


H 


0.7 


0.1 


0.6 


318 


319 


203 


K 


0.6 


0.1 


0.6 


296 


283 


208 


L 


0.4 


0.1 


0.4 


206 


179 


127 


M 


0.4 


0.1 


0.4 


227 


178 


162 


N 


0.4 


0.1 


0.3 


105 


103 


69 



Seyfert II nucleus, but MCG +0-29-23, NGC6810, and 
MRK897 do not. Contamination by a starburst com- 
ponent in NGC 34 and NGC 5953 seems likely. We will 
present data for all objects in this work, but exclude the 
three canonical HII nuclei from analysis in paper II. 

2.2. Imaging and data reduction 

The following galaxies in our sample have exist- 
ing high quality SEDs in lAlonso-Herrero et alJ (|2003f ): 
NGC 1068, NGC 4388, NGC 4968, NGC 5506, NGC 7172 
and NGC 7674. The remaining 42 Type II Seyfert galax- 
ies were observed in the following infrared bands: near- 
Infrared (NIR) J (1.25/im), H (1.65^m), and K (2.2/.tm) 
bands; and mid-IR (MIR) L {ZM^lm), M (4.66/im), 
and N (10/im) bands. Details of the observing runs are 
listed in Table |3l typical seeing conditions (for a Ho=70 
km/s/Mpc) are summarized in Table 01 

2.2.1. NIR observations 

Observations in the near-IR bands followed the usual 
strategy for this wavelength regime: several images of 
short integration time (3-20 seconds) are obtained, to 
avoid saturating the detector and to minimize fluctua- 
tions in sky brightness; the location of the object on the 
detector is changed for each image, to subtract the sky 
efficiently and to avoid bad pixels. 

The data were reduced, in part, with the XMOSAIC 
task from the IRAF XDIMSUM package. Darks were 
subtracted and flat fields applied. Sky images were cal- 
culated from a masked median of the object images or 
from the mean of sky-only images. 

Sky images were subtracted and results were registered 
and shifted before co-adding. The standard deviation 
per pixel was calculated from an exposure map. 

The observations with ISAAC in visiting mode were 
done under non-photometric conditions. Flux calibra- 
tion was achieved using field stars found in the 2MASS 
catalog^ as standard flux calibrators. This strategy did 

^ This publication makes use of data products from the Two 



not work for NGC 424, which had no field stars in the 
galaxy frame, and those observations were lost. 

2.2.2. MIR observations 

At these wavelengths the atmosphere is very bright, 
and temporarily and spatially variable; the telescope and 
surroundings also emit and refiect radiation. So obser- 
vations were made using standard chopping and nodding 
techniques. 

Mid-IR detectors are small, and the probability of find- 
ing a star in the same frame as the galaxy is slim, so each 
mid-IR galaxy observation was followed by observation 
of a star to determine the Point Spread Function (PSF). 
Known standards were used for flux calibration. 

The ISAAC mid-IR observations were reduced with 
the ECLIPSE "jitter" pipeline. In a number of cases ob- 
serving conditions degraded at the end of the integration, 
so some of the nodding series were discarded by hand. 

The Gemini N band data were reduced using the 
MIDIR task from the IRAF GEMINI package. When 
pipeline results were compared with a reduction by hand 
the images were very similar. 

3. CONSTRUCTING THE SEDS 

The SEDs contain nuclear magnitudes extracted from 
near- and mid-IR images. In t j3.1l we describe how 1-D 
brightness profiles are calculated from 2-D images; the 
modeling of these profiles, to separate the nuclear com- 
ponent, is described in tJ3.2l 

3.1. Determining the Surface Brigfitness Profiles 

Deep IR images of a sample of Type II Seyfert galax- 
ies were obtained. These images cover the entire galaxy, 
but only the nuclear emission is needed for the SED. 
To deconvolve the galaxy into its different components 
— bulge, disk, bar, and nucleus — we first deter- 
mine the ra dial Surface Brightnes s Profile (SBP) of 
the galaxy (Kotilainen et al. 1992'; Zitelh et al. '1993t 
lAlonso-Hcrrcro ct al. 1996; Prieto ct al. 2001, 2002). 
The main advantage of this 1-D approach is that smaller 
structures (spiral arms, HII Regions, rings) are diluted 
relative to the more fundamental structures (bulge, disk, 
nucleus and bar) and the quantity of data that must be 
fitted is reduced from ^500x500 pixels to ^200 points. 

The task ELLIPSE from the IRAF STSDAS package is 
used on the sky subtracted images to obtain the SBPs. 
A manual masking removes contamination from differ- 
ent sources that might distort the isophote fitting, like 
foreground bright stars or other galaxies. ELLIPSE fits 
ellipses to the galaxy isophotes; these are described by 
the center coordinates, the position angle and the ellip- 
ticity. The center coordinates were held fixed while the 
other parameters were allowed to vary. A linear radial 
sampling was used, with steps of 0.5-1 pixel, depending 
on the size of the galaxy. In cases where the edge of 
the galaxy was too faint and the ellipse fitting process 
did not converge, the fitted isophotes observed in other 
images were adopted. In order to take into account the 

Micron All Sky Survey, which is a joint project of the Univer- 
sity of Massachusetts and the Infrared Processing and Analysis 
Center / California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 
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inclination of the galaxy, the SBP was corrected using 
the mean ellipticity of the last fitted isophotes. 

3.2. Modeling the Surface Brightness Profiles 

Most spiral galaxies can be described using a disk, 
bulge, optional bar, and, sometimes, an active nucleus. 
Historically each component has been modeled with a 
specific analytical function. These are described below. 

Nucleus — The nuclear point source is modeled using 
Dirac's Delta: 



S{t)dt = l 



(1) 



So the nuclear SBP is described by: 

S„(r)=a„J(r), (2) 

where cr„ is the nuclear surface brightness. 

Disk — The disk is modeled using an exponential law: 

=CTdexp{-r/rd}, (3) 

where ad is the surface brightness of the disk at r = 0, 
(Trf = Srf(r = 0), and is the disk scale length. 

Bulge — The bulge is modeled using a general Sersic's 
law: 

SB(r) = (TBexp{-6„-(r/rs)^/"}, (4) 

where n is a free parameter and 6„ is defined by: 

r(2n)=27(2n,fe„), (5) 

where F is the (complete) gamma function, and 7(2n, x) 
is the incomplete gamma function defined by 



7(2n,x): 



dt. 



(6) 



Common values of 6„ are &i — 1.678, corresponding 
to an exponential SBP, and 64 — 7.6692, correspond- 
ing to a de Vaucouleur's profile. The most usual ana- 
lytical expressions for 6„ are 6„ = 1.9992 • n — 0.3271 
for 0.5 < n < 10, and &„ = 2n - 1/ 3 for n > 8 
(jGrahaml I2OOII: iGraham fc Driveii l2005f) . In this work 
we adopt the first approximation, because the typical 
values of n are bet w een 0.5 and 4 in the IR bands 
jBalceUs et al.' '2003; 'Hunt et al.' '2 0041 : lGrosb0l et all 
[2004; Kormcndy & Kennicutt 2004). 

Bar — Many galaxies in the sampled are classified 
as barred, and approximately 2/3 of all spiral galax- 
ies classified as non-barred look barred in the IR 
(|Kormendy &: KennicuttI (|2004l ). and references therein). 
In order to model the bar, a second general Sersic's law 
was included in the model, which introduced 3 new free 
parameters. 

Fmal Model — The final model is the sum of the 3 or 4 
components convolved with the image PSF in order to 
introduce the effect of the seeing: 



E(r) = [a-„(5(r) + a^e ''^ + ase 



PSF{r) (7) 




5 10 15 20 

Radius [arbitraiy units] 

Fig. 2.— Example of the Moffat fit of a PSF star. The observed 
data are the blaclc diamonds and the Moffat profile is represented 
by the red line. A Gaussian profile is shown for comparison (blue). 



3.2.1. The Point Spread Function 

To determine the correct parameters for the galaxy 
profile we need a precise representation of the PSF de- 
scribing the seeing and instrument response of each ob- 
servation. For our NIR imaging, the PSF is very accu- 
rately determined from a star found in the galaxy frame, 
which guarantees the same distortion for galaxy and star. 
In the MIR imaging the star is observed before and/or 
after the galaxy and, as the integration is quite long at 
these wavelengths, seeing conditions can change signifi- 
cantly. This means that the MIR stellar PSF may not 
accurately match the galaxy. In such cases the bulge 
profile is taken from the NIR fit, but scaled to match the 
external part of the observed profiles, and the nuclear 
flux was obtained assuming Fnuc = Fgai — Fbui ■ See Ap- 
pendix [X] for details. 

The PSF radial profile was obtained using the same 
technique as the galaxy profile (Section 13. 1|) and then 
modeled with a Moffat function to allow for interpola- 
tion / extrapolation when fitting the galaxy profile. The 
Moffat function is similar to a Gaussian, but has more 
extended tails which better fit the effects of seeing, as 
shown in Fig. [2l It is given by: 



Moffir): 



peak 



(l + (r/a)2)/3 



(8) 



3.2.2. The Fitting Routine 



To measure the nuclear flux — i.e., to deconvolve the 
nuclear emission from that of the galaxy — the radial 
SBP was modeled using the expression described above 
(Eqn. [T]) . The fitting routine uses the "Downhill Simplex 
Method in Multidimensions" first proposed by Nelder & 
Mead (1965): the minimization process starts with an 
initial guess from the user and then takes a series of steps 
exploring the available parameter space of the problem 
until it achieves the termination criterium which is, typi- 
cally, an insignificant variation between one step and the 
next. 

As the algorithm does not guarantee that the minimum 
found is global, it was started from 4 or 5 extreme points 
of the parameter space: the first initial guess includes all 
the galactic components of the model with comparable 
contributions to the total flux; other initial points each 
zero one of the components. These first guesses are not 
completely blind, but instead are based on the morphol- 
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ogy seen in the images. An iterative process that ran- 
domly perturbs the initial point is also included. Tests 
showed that 200 iterations were sufficient to reliably find 
the global minimum, which corresponds to the model 
that minimizes x^S the percentage contributions of the 
galactic components are obtained from this model. 

3.2.3. Errors in the fitted components 

The flux of each fitted galactic component has two 
main sources of uncertainty: the first from the fiux cali- 
bration (the zero points, obtained by the standard cali- 
bration procedure) and the second from the fitting pro- 
cess. The confidence limits for the fitted parameters, de- 
termined using the fitting routine described previously, 
are calculated using Monte Carlo simulations: the best- 
fit parameters are perturbed (one at a time) to con- 
struct many synthetic data sets. Each of these syn- 
thetic data sets has a particular Xsyn value and the first 
confidence level (corresponding to la error), which oc- 
curs 68.3% of the times, corresponds to the cases where 
Ax^ = Xmm ~ Xsyn — 1- The final errors are calculated 
as the quadratic sum of the individual errors for the two 
sources of uncertainty. 

3.2.4. Synthetic Data 

Synthetic data were used to test the fitting routines. 
SBPs of different "kinds" of galaxies were constructed: 
including all components equally; all components with 
very small or very large contributions; omitting one or 
more components; with wrong PSFs; and with different 
sampling steps (for the galaxy and the PSF star). 

The first test recovered the parameters of the synthetic 
SBP. The routine gives excellent estimates of the parame- 
ters for galaxies composed of 3 or 4 structural parameters 
with non-negligible contributions to the total emission. 
In cases when one or more components were not present, 
the results were also good, corresponding to observed ex- 
amples like the non-barred galaxies, or a couple of cases 
where nor disk nor bar is observed. Including all the 4 
components might give a better x^i because there will 
be more free parameters, but if the image does not show 
a bar, for example, it is obvious that it should not be 
included in the fit. 

Next, changes in the PSF were used to test the effects 
on the results (i.e., using a wrong PSF to describe the 
seeing conditions). We found that the results could be 
recovered (this is, at least 2 of the initial starting points 
could reach the actual solution) whenever the FWHM 
was narrower or broader by up to 50% of its real value. 
These tests apply especially to the mid-IR observations: 
as the integration time is long, the probability of a change 
in the seeing conditions between the observation of the 
PSF star and galaxy is greater. 

The truncation of the galaxy profile severely affects 
the performance of the fitting routine. Fortunately there 
are very few such cases in the observations, and then, 
typically, only one of the broad-band images is affected, 
which means that results can be checked with the unaf- 
fected profiles and parameters constrained appropriately. 

Finally, it is found that the sampling step is not par- 
ticularly important, providing that the PSF is sampled 
by a finer discretization than the sampling of the galaxy. 

3.3. Messier 31 



TABLE 5 

Values for the relevant fitted parameters of M 31 
decon volution process. 



Component 


1.25/jm 


1.65/im 


2.2^tm 


3.6^m 


4.5^m 


5.8/jm 


Bulge r^ff 
Bulge n 


185 
2.0 


228 
2.3 


207 
2.2 


301 
2.5 


259 
2.4 


209 
2.2 


Disk r 








1165 


991 


996 


Bar r^ff 
Bar n 








1.2 
0.6 


1.2 
0.5 


0.6 
1.2 



Note. — The scale lengths are expressed in arc-seconds. The bright 
excess in the central part of the galaxy profiles is modeled as a bar in 
the IRAC bands and as a nuclear component in the 2MASS bands 
(due to the differences in pixel scale). The disk is not included in the 
2MASS fits because of the smaller image coverage. See text for details. 



To test the SBP fitting procedure. Messier 31, the An- 
dromeda galaxy, was studied in the same way as the ob- 
jects in our sample. 

Images of M31 were taken by IRAC as part of the 
Spitze r General Observer p rogram, ID 3126, in August 
2005. iBarmby et all ([20061 ) measured the SBPs by fit- 
ting elliptical isophotes to the images using the ELLIPSE 
task. Profiles in the three longer wavelength bands (at 
4.5, 5.8, and 8.0 fim) were measured using the ellipses fit- 
ted to the 3.6 fim image, so that colors could be measured 
at the same spatial locations from all 4 ban ds. These pro- 
files w ere kindly provided by P. Barmby (IBarmbv et al.l 
l2006h . 

Images in the J, H, and K bands were obtained from 
the 2MASS public catalog. The profiles were constructed 
using ELLIPSE. However, M31 is so close that the 
2MASS images do not cover the full extent of the galaxy: 
while Spitzer-IRAC covers > 4000" in radius, the images 
of 2MASS cover only ~500". This compromises the re- 
sults for the disk component in the 2MASS data. 

M31 is one of the nearest galaxies, and all galactic 
components can be detected in all the IR bands, which 
is not the case for the galaxies in our sample. The SBPs 
of M31 were fitted from 1.25/im to 5.8^j.in (Fig. E]). The 
8.0/im SBP was not included because it is dominated by a 
10 Kpc star-forming ring and individual stellar features. 
The 5.8/im SBP was truncated at ^1300" to avoid wig- 
gles produced by the 10 Kpc ring. The relevant fitted 
parameters are listed in Table [5j 

It is interesting to note the bright excess in the cen- 
tral part of the galaxy profiles. This excess is modeled 
as a bar in the IRAC bands and as a nucleus i n the 
2MASS bands (due to the difference in pixel scale). iPengj 
([2002) studied the central region of M31 in detail us- 
ing the /bubble Space Telescope (HST). He modeled the 
center of the galaxy using 6 components to account for 
3 different peaks of emission in UV and optical images. 
However the UV peak disappears in the near-IR images 
and the 2 V-band peaks diminish, so here we assume that 
the nucleus/bar fitted by the process is real, but the pixel 
scale does not allow a more detailed modeling. 

As the disk could not be fitted in the 2MASS profiles 
only the fitted radius of the bulge is compared: 700 pc 
in the J band, 862 pc in the H band, 787 pc in the K 
band, 1140 pc at 3.6 fim, 983 pc at 4.5 fim, and 795 
pc at 5.8 fim. The dispersion of these values about the 
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Fig. 3. — Radial surface brightness profile fits of J, H, and K bands (top panels), and 3.6/jm, 4.5^m, and 5.8/xm (middle and lower 
panels) of galaxy Messier 31. The abscissas show Log(counts/pix^) and the ordinates show the distance from the center in arc-seconds. 
The percentage residues are shown in the bottom panels. The model (red thick continuous line) is composed by a nuclear PSF (pink 
dotted-line), a bulge (green dashed-dotted line), a disk (blue dashed line), and a bar (yellow dashed-double-dotted). The black points 
correspond to the observation. The brightness excess in the central part of the galaxy is modeled as a bar in the IRAC bands and as a 
nuclear component in the 2MASS bands. See text for details. 



mean is between 20 and 30% and they show no significant 
systematic change moving toward longer wavelengths. 

The results corroborate some naive expectations: un- 
less there are strong color gradients within each of the 
structural components (the bulge, disk, and/or bar) the 
scale lengths should remain fairly constant for the galax- 
ies in our sample, while the brightness scale factor should 
decrease from the NIR to the MIR as the stellar popu- 
lations become dimmer. Also, since the bulge has older 
stars and is redder than the disk and the bar, it is ex- 
pected to be the brightest component in the mid-IR. 

Given these findings for M31, we decided, when fit- 
ting the SBPs of the galaxies in our sample, to con- 
strain the scale lengths in the MIR using the values fit- 
ted in the NIR, since the very dim stellar components 
under-determine the fitting at these wavelengths. See 
Appendix [K\ for details. 

4. GALAXY FITTING RESULTS 

The fitting of the SBPs for each of the galaxies in our 
sample was done assuming the model described in §3.21 
Visual inspection of the images and the Hubble classi- 
fication of each galaxy were used to check for coherent 
results. 

Fitted parameters are generally consistent between ob- 
serving bands. The following is a brief summary of the 
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Fig. 4. — Contribution from the nuclear flux to the total flux as 
a function of band: J, K, L and M measurements are presented 
with triangles, squares, pentagons, and circles, respectively. The 
right-hand side histograms show the projected distributions. 
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Fig. 5. — Spectral Energy Distributions determined for our sample. Note the different scale in Y axis i n galaxies NGC 1068, MCG -2-8-39, 
NGC 1194, IRAS 03362-1642, NGC 5506, and NGC 7590. 6 SEDs are from I Alonso-Herrero et~all ((2003). 



results, which are also presented in Table IHl 

No disk or bar was detected in the M band, while the 
bulge was detected in 7/28 sources. In those objects the 
bulge contributed a mean brightness of 44% to the total 
flux. Only one bar was detected in the L band with a 
contribution of 17%. In the NIR bulge, disk and bar (if 
present) have similar contributions of ^ 30-50%. 

Seven out of 105 near~IR fits require a null nuclear 
contribution. The mean percentage contribution of each 
galactic component is summarized in Table [H A mean 
nuclear contribution of < 10% is found in the NIR bands 
(JHK), while a strong increase can be seen in the MIR, 
going from 40% in the M band to - 100% in the M 



and N bands. 

The percentage contribution of the nuclear component 
in the J, K, L and M bands is shown in Fig. 0] as a func- 
tion of the OIII luminosities (tabulated in Paper II). It 
can be seen that, as expected, the percentage nuclear 
contribution increases for more luminous AGN, as indi- 
cated by the OIII measurements (Kotilainen et al. 1992; 
Alonso-Herrero et al. 2001, 2003). 

The fitted parameters for the galactic structures are 
tabulated in Table [T] In Appendix |X] details of the in- 
dividual SEP fittings are given. The fits for the MIR 
profiles of all the galaxies are shown at 2 different scales: 
a detailed one (bottom panels), showing the central part 
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Fig. 5. — Continuation. 
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of the galaxy, and a global one (middle panels), where 
the whole galaxy is observable. 

4.1. Spectral Energy Distributions 

The constructed nuclear IR SEDs are shown in Fig. [Sj 
The nuclear fluxes and standard deviations for each ob- 
served band are summarized in Ta ble El Here we include 
the SEDs already determined by lAlonso-Herrero et al.l 
(|2001, 200S)- Upper limits when the galaxy was unde- 
tected in L, M, and N bands are also shown. 

There are some cases where a lack of photometric 
points prevents us from obtaining a good SED, as with 
NGC 1194, IRAS 03362-1642, ES0 33-G2, NGC 4941, 
and NGC 7590. The SEDs of both nuclei in IRAS 00198- 
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7926 are similar and do not determine whether the north- 
ern or southern nucleus (or both) corresponds to the 
AGN. 

Visual inspection of the SEDs reveals that we can 
broadly group most of them into two types depending 
on the NIR behavior: one group shows a "classic" bell- 
shaped SED, as expected from the emission of a dusty 
torus, while the remaining SEDs show a much flatter 
NIR spectral distribution, with a clear excess with re- 
spect to the first group. We have quantified this behav- 
ior by measuring the NIR slope of the SEDs and classify 
them depending on whether the slope a (A/a oc A") is 
larger or smaller than 1. Individual and average SEDs for 
these two groups can be seen in Figure [6l The tabulated 
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TABLE 6 

Mean, standard deviation, skewness, median and mode 
percentage contribution of each galactic component 



Band Number Mean a Skew Median Mode 
(%) (%) (%) (%) (%) 



100 
44 
13 
0.8 
0.8 

44 
57 
35 
33 
34 



36 
47 
52 
55 



17 
28 
30 
30 



Note. — In column 'Number' we indicate the number of fits which 
include that particular component over the total number of images 
available in a particular band. Results obtained by AH03 are not 
included in this Table. 

values for the NIR slope are shown in Table \8\ Notice 
that some values have rather large errors and therefore 
the classification can be ambiguous. The possible nature 
of the SEDs presenting a NIR excess is discussed in the 
next section and further analysis can be found in Paper 
II. 

In the same figure we also include the SEDs of the 
three HII nuclei found in our sample and of NGC 34 and 
NGC 5953 where contamination by a starburst compo- 
nent is likely. It can be seen that the morphology of their 
SEDs clearly departs from the bell-shaped spectra of the 
first group of Seyfert nuclei and are somewhat closer to 
the second group of Seyfert galaxies, as it is also noted 
by the flat value of their NIR slopes (see Table [8]) . 

4.2. Near-IR excess 

In a Type I Seyfert, where the inner region of the cen- 
tral source is directly observed, the accretion disk can 
make a significant contribution to the total near-IR emis- 
sion. Recently Landt et al. (2011) and Lira et al. (2011) 
have shown that this is necessary to explain the contin- 
uum around 1/im in these objects. 

A "typical" Type II source, on the other hand, should 
have an SED that decreases monotonically towards 
shorter wavelengths. But inspection of the SEDs in this 
work shows several NIR SEDs with a clear excess which 
resembles Type I systems (see Fig. [5] and Table |9]). 

All galaxies in our sample had an early classifica- 
tion as Type II systems (with MCG -3-34-64 being clas- 
sified as Seyfert 1.8, NGC 5506 as a Seyfert 1.9, and 
NGC 7130 as a LINER). However, as noted above, more 
recent observations show that NGC 6810, MRK 897, and 
MCG -1-0-29-23 correspond to HII nuclei, while NGC 34 
and NGC 5953 might have contamination by a starburst 



component. An IR excess is expected for these ob- 
jects because, in the absence of an active nucleus, the 
Rayleigh- Jeans emission tail from cold stars dominates 
this spectral region for moderately old starbursts (~ 
a few dozen Myrs). This is clearly observed in other 
starbursts and extensively modeleded (e.g., Efstathiou 
et al. 2000; Dopita et al. 2005; Rodrfguez-Merino et 
al. 2011). 

The remaining Seyfert II objects that show an IR ex- 
cess require alternative explanations; we discuss various 
mechanisms below. 

A traditional torus — a spatially continuous struc- 
ture with large optical depth — would completely ob- 
scure the accretion disk in Type II objects, but a clumpy 
torus gives a non-zero probability of observing the central 
source at any orientation; the latter scenario is increas- 
ingly accepted and could explain some near-IR emission. 

Alternatively, a nuclear jet could make a significant 
contribution to the observed near-IR fluxes. In this case 
the emission would be due to synchrotron instead of ther- 
mal processes. Jet emission is found to be prominent 
in radio-loud AGN and in relativistically boosted radio- 
quiet AGN (also known as radio-intermediate AGN; Fal- 
cke et al. 1996; Barvainis et al. 2005). Radio-loud AGN 
are defined as those with Rl > 100 (where Rl = 
FbGHz/ Fb\ Kellermann et al. 1989), radio- intermediate 
AGN are those with 3 — 10 < Rl ^ 100, while radio-quiet 
AGN have i?^ << 10. 

We have compiled 5 GHz measurements for our sam- 
ple and determined Rl by computing Li°°^ from mea- 
sured [OIII] fluxes and then turning them into nuclear B 
fluxes following Marconi et al. (2004). The resulting Rl 
values are presented here for the SEDs with NIR excess 
and in Paper II for the complete sample. We found that 
~ 30% of the sources with NIR excess can be classified 
as radio-loud (F 01475-0740, NGC 4501 and NGC 7496); 
only two additional radio-loud sources present a bell- 
shaped SED (F 04385-0828 and NGC 7172), representing 
a 10% of that subsample. NGC 4501 presents a rather pe- 
culiar SED, with a power-law distribution across the en- 
tire l-lO/im range (see Paper II). In this case synchrotron 
emission might explain the spectral shape. 

In summary, jet emission might explain some observed 
features in the infrared SEDs, but the link is not proven 
and does not explain all cases. 

A very hot dust component has been proposed to 
explain near-IR emission in Mrkl239 and Mrk766 
(Rodriguez- Ardila et al. 2005; Rodriguez- Ardila & Maz- 
zalay 2006). These are type I AGN, which provide an 
unimpaired view of the hottest dust region, and emission 
from this dust — possibly carbonaceous grains surviving 
inside the silicate sublimation radius — appears as an 
excess superimposed on the disk component. Interest- 
ingly, the peak observed in these two sources is clearly 
distinct from the putative torus component (which is not 
observed in these sources as no data are available beyond 
4/im). 

Similarly, a hot dust component was used in other 
works when modeling the near-IR region of Spitzer spec- 
tra for powerful QSOs (Mor & Trakhtenbrot 2011, Mor et 
al. 2009, Schweitzer et al. 2008; Deo et al. 2011). They ar- 
gue that this component is required to fit the NIR part of 
the spectrum and that disk emission can only contribute 
significantly at l/um, but not enough at longer wave- 
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TABLE 7 

Fitted parameters of the galaxy deconvolution in each band of observation. 



Galaxy 
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" nuc 

Watts/m J 


"Bui 
Watts/m J 


eff 

pc 




Watts/m //^/[ J 


^ dis 
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"bar 
Watts/m /;U/[ J 


eff 

pc 




NGC34 


J 


5.6e-14 
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Note. — 



The remaining of this Table is available in the online version of the Journal. 
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TABLE 8 
Nuclear Fluxes 



Galaxy 




J band 


H band 


K band 


L band 


M band 


N band 


NIR slope h 






mjy 


mJy 


mJy 


mJy 


mJy 


mJy 




NGC 34 


1 


.60±0.25 


5.17±0.28 


4.92 ±0.70 


1.27 ±0.17 


1.19 ±0.16 


171 ± 12 


0.97±1.26 


IRAS 00198-7926 N 





.35±0.02 


0.44±0.05 


0.65 ±0.08 


13.10 ±1.21 


9.27 ±0.75 


— 


0.10±0.15 


IRAS 00198-7926 S 





.llitO.Ol 


0.54±0.08 


3.28 ±0.43 


13.10 ±1.21 


9.27 ±0.75 


— 


5.02±0.16 


IRAS 00521-7054 


2 


.57±0.14 


4.70±0.29 


11.66 ±1.65 


19.91 ±1.84 


39.31 ±2.79 


- 


1.69±0.28 


ES0 541-IG12 


1 


.35±0.16 


4.56±0.56 


8.73 ±1.44 


19.15 ±3.18 


15.71 ±2.08 


<140 


2.28±0.62 


IRAS 01475-0740 





.66±0.02 


1.00±0.04 


1.41 ±0.19 


3.57 ±0.45 


8.47 ±0.63 


- 


0.34±0.09 


NGC 1068t 




9.8 


97.6 


449.6 


3691 


8245 


23000 


5.77±0.86 


NGC 1144 





.35±0.11 


0.34±0.15 


0.42 ±0.14 


1.84 ±0.04 


5.10 ±0.55 


23±10 * 


-0.67±0.24 


MCG -2-8-39 





.07±0.05 


0.07±0.09 


0.07 ±0.23 


1.25 ±0.19 


5.39 ±0.47 


<114 * 


-1.00±0.01 


NGC 1194 


0, 


.02±0.12 


0.17±0.23 


1.01 ±0.90 


- 


- 


133 ± 19 


5.92±0.43 


NGC 1320 





.llitO.lO 


0.35±0.65 


0.98 ±0.97 


- 


- 


195 ± 27 


2.86±0.17 


IRAS 03362-1642 





.04±0.06 


no-det 


no-det 


- 


- 


970 ±110 


- 


IRAS 04385-0828 





.19±0.10 


0.90±0.14 


2.36 ±0.84 


- 


- 


254 ± 46 


3.43±0.64 


ES033-G2 





.36±0.38 


1.47±0.57 


3.21 ±0.80 


- 


- 


- 


2.86±0.68 


IRAS 05189-2524 


1 


.93±1.18 


7.21±3.37 


14.67 ±5.23 


- 


- 


463 ± 62 


2.57±0.66 


MCG -1-0-29-23 





.42±0.13 


0.53±0.22 


0.75 ±0.04 


1.33 ±0.22 


10.13 ±1.22 


73 ± 7 


0.02±0.11 


NGC 3660 





.74±0.05 


1.07±0.21 


1.09 ±0.23 


1.26 ±0.07 


<11.07 


24.2±3.3 


-0.32±0.37 


NGC 438St 




0.06 


0.71 




39.9 




245±15 


4.18±2.17 


NGC 4501t 


4 


.58±1.12 


4.93±1.18 


4.76 ±0.71 


2.68 ±0.25 


<7.22 


6.0±0.5 * 


-0.93±0.11 


TOL 1238-364 


1 


.67±0.30 


1.96±0.27 


2.27 ±0.64 


4.81 ±0.55 


7.96 ±0.68 


31.2±0.4 


-0.46±0.02 


NGC 4941 




no-dct 


no-det 


no-det 




<12.23 


47 ± 5 




NGC 4968t 




no-dct 


0.6 


3.7 


23.2 


56.5 


280±20 


5.32±0.01 


MCG -3-34-64 





.81±0.14 


1.13±0.54 


1.32 ±0.34 


8.29 ±0.34 


13.75 ±0.96 


603±30 


-0.14±0.19 


NGC 5135 





.17±0.45 


0.19±0.38 


1.74 ±1.41 


3.88 ±0.10 


8.18 ±0.62 


104 ± 10 


3.16±2.08 


NGC 5506t 




13.8 


59.0 


120.4 


340.1 


530.0 


100 ± 10 


2.83±0.80 


NGC 5953t 


1 


.61±0.53 


2.77±0.18 


1.06 ±0.05 


0.23 ±0.02 


<12.23 


24±7 


-1.68±1.53 


MCG -2-40-4 


3 


.57±1.47 


7.10±3.00 


15.99 ±6.15 


32.63 ±2.21 


21.43 ±2.24 


355± 21 


1.65±0.10 


IRAS 15480-0344 





.96±0.09 


1.62±0.62 


2.24 ±0.18 


4.95 ±0.17 


19.11 ±1.68 


100 ± 10 


0.50±0.23 


NGC 6810 


1 


.31±1.70 


1.30±8.41 


1.06 ±0.84 


1.10 ±0.18 


6.23 ±0.52 


49± 6.5 


-1.37±0.21 


NGC 6890 





.17±0.22 


0.24±0.42 


0.86 ±0.73 


4.16 ±0.27 


8.61 ±0.48 


107± 14 


1.89±0.92 


IC 5063 





.82±0.95 


0.14±0.79 


1.20 ±0.50 


16.38 ±1.17 


52.48 ±4.39 


787±105 


-0.37±4.00 


MRK897 





.43±0.16 


0.75±0.45 


0.79 ±0.11 




<11.15 


5.7±0.8 


0.08±0.54 


NGC 7130 


1 


.53±0.50 


2.32±1.56 


2.58 ±1.15 


4.41 ±0.41 


9.27 ±0.79 


197± 26 


-0.08±0.33 


NGC 7172t 


non-dct 


<0.4 


3.4 


30.0 


61.4 


103 ± 10 




MCG -3-58-7 


1 


.45±0.59 


4.47±2.31 


10.12 ±4.11 


20.29 ±1.07 


41.19 ±4.94 


218 ± 15 


2.42±0.35 


NGC 7496 


1 


.07±0.40 


1.77±1.44 


2.36 ±3.31 


3.16 ±0.13 


21.28 ±3.55 




0.39±0.24 


NGC 75820 





.69±4.33 


1.64±21.24 


5.45 ±19.84 


79.06 ±3.88 


63.19 ±4.72 


195 


2.66±0.31 


NGC 7590 




no-det 


no-det 


0.26 ±0.55 


no-det 








NGC 7674t 




1.25 


5.0 


12.3 


53 




344 


3.04±0.54 


CGCG 381-051 





.91±0.19 


1.01±0.34 


1.37 ±0.43 


<1.41 


<11.65 


75 ± 8 


-0.27±0.20 



Note. — f indicates when most of the SED measurements were determined by AHOl and AH03 (see Table 1 for 
more details); * indicates data from Joy & Ghigo (1988) and Maiolino et al. (1995). Caution must be taken because of 
the very different dates of the observations when compared with our shorter-wavelength data. J N-band observations 
from Gemini and G04 only provided upper limits, t] NIR slope a {Xf\ oc A°) was tabulated using least square fits to the 
data. Photometry for NGC 7582 will be complemented with observations from Prieto et al. 2010 and Ramos- Almeida 
et al. 2009 in paper II. 
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Fig. 6. — Normalized Spectral Energy Distributions for three groups of sources. The left-hand-side panel shows nuclei with a NIR slope 
> 1. The central panel shows those sources where a NIR excess is present, i.e., with a NIR slope < 1. Finally, the right- hand- side panel 
shows the SEDs of the three HII nuclei found in our sample and two nuclei which are likely to be contaminated by nearby starburst 
emission. Individual and average SEDs are shown for the first two cases, ignoring upper limits when deriving the averages. 



TABLE 9 
Near-IR Excess and Radio-loudness 



Galaxy 


Excess 


HII Nuclei 


Rl 


IRAS 00198- 7926 N 


yes 


no 




IRAS 01475-0740 


yes 


no 


70 yes 


NGC 1144 


yes 


no 


3.4 no 


MCG-2-8-39 


yes 


no 


0.2 no 


NGC 3660 


yes 


no 


0.9 no 


NGC 4501 


yes 


no 


95 yes 


Tol 1238-364 


yes 


no 


0.4 no 


MCG-3-34-64 


yes 


no 


2.2 no 


IRAS 15480-0344 


yes 


no 


0.6 no 


NGC 7130 


yes 


no 


1.1 no 


NGC 7496 


yes 


no 


32 yes 


CGGC 381-051 


yes 


no 


2.4 no 



lengths. Therefore very hot pure-graphite dust emission 
is proposed. 

A similar very hot dust component could explain some 
of our observations, but, because of its location inside 
the torus sublimation radius, it again requires a clumpy 
absorbing medium to explain its detection. 

A very compact nuclear starburst might have a signif- 
icant contribution to the nuclear fluxes. To test this hy- 
pothesis we have looked at the positions of Seyfert nuclei 
in the diagnostic diagrams presented in Figure [TJ where 
nuclei with a NIR excess are shown with a square. It can 
be seen that there is no clear segregation of these nuclei 
towards the HII location in the plot. Also, these nuclei 
do not show preference for the presence of a very young 
stellar continuum as shown in Table [2j 

Finally, a young and luminous nuclear stellar cluster 
could contribute to the observed near-IR emission. How- 
ever, with an absolute magnitude zab = ~13 for the 
most luminous examples (Cote et al. 2006) their flux is 
not sufficient to make a significant contribution to the 
observed SEDs. 

5. SUMMARY 



We have constructed nuclear IR SEDs of 40 Type II 
Seyfert galaxies from ground-based, high-resolution ob- 
servations. The images were collected between 2003 and 
2004, to reduce variability distortion of the SEDs as much 
as possible, and under good seeing conditions. 

We also obtained optical spectroscopy for most of 
the sample. The spectra were modeled using the 
STARLIGHT code to characterize their stellar contin- 
uum properties and emission lines were used to classify 
the spectra as active or HII nuclei. Three objects from 
our sample most likely harbor a starburst nucleus and 
not an AGN, in agreement with other findings. 

The surface brightness profiles of each galaxy, in each 
observed band, were constructed by fitting elliptical 
isophotes to the images. The profiles were then mod- 
eled using a nucleus, bulge, disk and, where necessary, a 
bar. This procedure was tested using synthetic data and 
IR observations of the Andromeda galaxy, which showed 
that the scale length of the different galactic structures 
does not change when observing at longer wavelengths. 
This allowed us to fix the scale lengths in the MIR 
when analyzing undetected components of distant galax- 
ies. Careful modeling of the PSF was also included to 
take into account seeing distortions. 

The fitted SEP parameters for any particular object 
are generally consistent between bands. In some cases, 
as in IRAS 03362-1642, NGC 4941, and NGC 7590, no 
nuclear component was detected by the deconvolution 
process in the NIR. 

A variety of SED shapes is found in our sample, despite 
the common Type II Seyfert classification. About 40% of 
the Seyfert SEDs are characterized by an upturn or ex- 
cess in the near-IR. The 5 objects with strong starburst 
emission (3 HII nuclei and 2 Seyfert nuclei with strong 
contamination from nearby starburst regions) also show 
a similar excess. For genuine AGN this component could 
be explained as emission from a jet, the accretion disk, 
or from a very hot dust component, leaking out from 
the central region through a clumpy obscuring structure. 
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The presence of a very compact nuclear starburst as the 
origin for this NIR excess emission is not favored by our 
data. 



We are forever indebted to Roberto Cid Fernandes for 
his help using the STARLIGHT code. LV gratefully 
acknowledges fellowship support by project MECESUP 
UCH0118, and partial support from Fondecyt project 
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APPENDIX 

NOTES ON THE SURFACE BRIGHTNESS PROFILE DECONVOLUTION OF INDIVIDUAL SOURCES 

In the notes below, whenever a parameter is described as "constrained" by some other value, it means that the 
parameter was allowed to vary within a limited range around the constraining value (s). Without this approach the 
parameter would adopt a solution that was non-physical or at odds with values found at other wavelengths. 

NGC34 — In this object it is clear how the disk becomes weaker compared to the bulge and nucleus from the J to L 
bands. The disk in the L band had to be constrained to have a similar scale length to the J, H and K bands. In the 
M-band the disk is not visible, but the bulge is still dominant. Here the bulge effective radius had to be restricted to 
have similar values to those derived at other wavelengths. The fits are shown in Fig. \7\ 

IRAS 00198-7926 — The NIR images for this galaxy show two very close compact objects. It is therefore very difficult 
to determine which of the two corresponds to the active nucleus. We label them as N and S. Elliptical photometry was 
performed on both objects in each band. The fitting results are shown in Fig. [7] (online). The southern object has a 
strange behavior, where the effective radius and the exponent of the Sersic's law of the bulge continuously decreases 
from J to K. Only one object is detected in the L band image, and it is included in both SEDs (north and south) 
because it is impossible to determine its correspondence with any of the two nuclei. Also, the L-band observations 
seem to suffer from seeing variations, as can be inferred from the fitting results: constraining the parameters of the 
bulge to be similar to the NIR values produces a very poor fit. As a result, all the emission observed in this band was 
considered as coming from the nucleus. 

IRAS 00521-7054 — This galaxy is classified as an E-SO but, in order to obtain a good fit, the disk component had to 
be excluded. Because the seeing changed between the observation of the PSF star and the galaxy in the L band only 
the external parts of the observed profile were considered when fitting the bulge (using the values found in the NIR 
bands) and the nuclear fiux was obtained assuming Fnuc = Fgai — Fbui ■ 

ESO 54I-IGI2 — This is a "text-book" active galaxy, with a bright dominant nuclear source. However, in the L-band 
the Sersic's index was constrained by the values obtained in the NIR fits, and in the M-band the effective radius of 
the bulge was constrained to the values obtained in the NIR and L-band. 

NGC424 — The NIR observations could not be flux calibrated and the MIR observations have no PSF star. Hence, 
no SED was obtained. 

IRAS OI475-O74O — In order to obtain a reasonable fit to the NIR profiles the scale length of the disk was constrained 
to be larger than that of the effective radius of the bulge. In the L band, the size and index of the bulge were restricted 
to have similar values to those found in the NIR fits. 

NGC1125 — This galaxy was only observed in the L and M bands. The disk and bulge parameters could not be 
restricted and the parameter degeneration was too large. Hence, no SED was obtained. 



18 



Videla et al. 




2 4 6 8 10 12 14 16 U 
Radius (arcsec) 
1 




4 6 8 10 12 14 16 1 
Radius (arcsec) 



4 6 8 10 12 14 16 18 
Radius (arcsec) 



L band, detail 





2 3 4 5 
Radius (arcsec) 




Radius (arcsec) 



Fig. 7. — Radial surface brightness profile of the J, H, K, L and M band observations of NGC34 (black dots with error bars). The y-axis 
corresponds to Log(counts/pix^) and the ordinates show the distance from the center in arc-seconds. The best-fit model is shown with a 
thick continuous line. Residuals in percentage are shown in the bottom part of each panel. The different physical components of the model 
are: the nucleus (dotted-line), the bulge (dashed-dotted line), the disk (dashed line) and a bar whenever necessary (dashed-double-dotted 
line). The mid-IR profiles (L, M and N, if available) are plotted in two scales: a detailed scale showing the central part of the galaxy 
(bottom panels), and a global scale, showing the complete galaxy (middle panels). Color plots for this and the remaining galaxies are 
available in the online version of the Journal. 
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Fig. 8. — Left: Raw image in the K band of galaxy IRAS 05189-2524. The obscured central pixels are evident. Right: Comparison of 
PSF star (dotted-blue Hne) and galaxy (dashed-red line) radial surface brightness profiles of L band of galaxy NGC 7590. It is evident that 
the galaxy is extended unless a strong variation in seeing conditions degraded the galaxy image to a level where no point source can be 
detected. Adding the fitting results in the near-IR and the resulting Spitzer spectra after star-formation subtraction, it seems this galaxy 
has at most a very weak AGN. 



NGC 1144 — This is an interacting system (NGC 1143/44, Arp 118) and the disk of the Seyfert galaxy is very distorted. 
Comparing the J, H, and K band images it is evident that the K band is less contaminated by the tidal tails, yielding 
well behaved parameters for the disk. Using these values to constrain the disk parameters in the J and H bands fitting 
results are consistent in all 3 bands. Results are plotted in Fig. [7] (online) and clearly show excess emission in the 
outer parts of the galaxy in the J and H bands, coming from the streaming material. In the L band all 3 components 
are discernible and the parameters are consistent with those of the NIR once the disk size is constrained to be larger 
than the bulge size. The galaxy is point-like in the M band, so all the emission is considered as nuclear and no fitting 
is needed. 

MCG -2-8-39 — This is a very nice example of a barred galaxy. The NIR fits were obtained without any restriction, 
and the results are consistent. The L band profile was fitted without disk or bar due to the bulge dominating the 
central region; the size and index of the bulge had to be constrained using the NIR values. 

NGC 1194 — This galaxy has a highly concentrated bulge [ns ~ 3.3). The results were obtained constraining the 
index value of the K band profile to that found in the J and H band fits. 

NGC 1320 — This galaxy was fitted limiting the effective radius of the bulge in the K band. 

IRAS 03362-1642 — The resulting fit of this galaxy does not include a nuclear component in the H and K bands. In 
the J band a weak nucleus appears, with a contribution of < 1%, almost negligible. The bar is very dominant in all 3 
bands. 

IRAS 04385-0828 — To fit this galaxy it was necessary to constrain the bulge Sersic's index in the J and H bands to 
the value found in K. 

NGC 1667 — This galaxy has a very complex morphology, which prevents the determination of a reliable profile fitting 
in the 3 bands observed. Hence no SED was obtained. 

ES0 33-G2 — No restriction was needed in J and H bands to fit this galaxy. In the K band the disk scale length was 
constrained to be larger than the bulge and bar sizes, obtaining consistent results for all the parameters. 

IRAS 05189-2524 — In order to obtain a consistent fit to this galaxy the disk size in the H band was constrained to be 
larger than the bulge size. In the K band the first point of the radial profile is lower than the second and third points. 
This is due to obscuration that can be observed in all the raw images in this band, as shown in the left panel of Fig. [Sj 
this first point is not included in the fitting process. In the N-band profile all the emission comes from the nucleus. 

ESQ 253-G3 — This is an interacting system whose two nuclei are too close to fit reliably. Hence, no SED was obtained. 

MCG -1-0-29-23 — Fig. [7] (online) shows the fitting results for this very special object. The results barely change from 
one band to another and the galaxy is very extended even in the L band. It is not clear whether the central spheroidal 
component corresponds to a bulge, since the images and the fitted Sersic's index both resembles more closely a bar. 
In any case, bulge and/or bar, the contribution to the central emission can be determined, even if it is not possible to 
unveil the exact nature of the component. 



20 



Videla et al. 



NGC3660 — This galaxy is very large, and due to the small size of the detector the galaxy has been truncated in the 
J and H bands. In order to determine the real value of the disk parameters the K band image was used and then these 
fitted values were used to constrain the disk parameters in the J and H band fits, excluding the profile for r > 45". In 
the L band, the profile with r > 2" was skipped because it is mostly sky contribution. The parameters of the bulge 
were restricted to those obtained in the J, H and K fits. 

NGC4501 — This is a huge galaxy compared to the detector, so the PSF was taken from the standard star. The only 
restriction needed to obtain consistent results was that the disk size was larger than the bulge size in the J and H 
bands. In the L band the change in seeing conditions is evident, and to overcome this problem the bulge was fitted 
for r> 3", constraining only the Sersic's index; the nucleus is then computed as Fnuc = F„ai — Fhui- The g alaxy was 
not detected in the M or N bands. The upper limit in the N band is that obtained by iGoriian et al.l ()2004[ ) 

TOL 1238-364 — The fit to this galaxy in the NIR bands does not include the bump produced by the ending of the 
flocculent arms between 16" and 23". In the J and H bands a disk scale length larger than the bulge size was imposed. 
In the K band the bulge index was constrained to have similar values to those found in the J and H bands. In the L 
and M bands the bulge and disk parameters were constrained to values found in the NIR fits. It seems that seeing 
conditions changed between the observations in the M band, since the PSF star did not produced an acceptable fit for 
the nucleus, so the bulge was fitted and the difference between the galaxy and the bulge was considered as the nuclear 
contribution, i.e., F^uc = Fgai - Fbui- 

NGC494I — This is a very big galaxy, which means that, again, the PSF is not contained in the galaxy frame but 
comes from observations of the standard star. Unfortunately the standard star is not a good approximation to the 
galaxy PSF, resulting in a null detection of a nucleus in the 3 NIR band fits. However, the galaxy was also not detected 
in the L and M bands, so it may also be that no detectable near- and MIR emission is associated with the active 
nucleus. 

MCG -3-34-64 — Without special constraints, the fitted NIR parameters are consistent. In the L and M band the 
seeing between the observation of the galaxy and the observation of the PSF star seems to vary, resulting in a poor fit 
in the L band. In order to fit both profiles, the bulge was fitted with the values found in the NIR fits and the nuclear 
contribution was calculated as Fnuc — Fgai — F^ui ■ 

NGC5135 — The bar in this galaxy is very strong and the disk is not seen in the IR images. Between I" and 3" there 
is a bright star-forming region, dominating the profile at all wavelengths; this region was excluded from the fit in all 
the profiles. The galaxy profile in the M band looks slightly broader than the PSF star, but the bulge component 
is not able to account for it, so a change in the seeing conditions is assumed to be responsible. The nuclear flux is 
calculated from the integration of the galaxy profile for r < I". 

MRK463 — This is an interacting system whose two nuclei are too close to fit reliably. Hence, no SED was obtained. 

NGC5953 — The bump observed near r=6" was excluded because it is produced by the ending of the flocculent arms. 
In the L band, the seeing is worse in the galaxy observation than the PSF star, as can be seen in Fig. [7] (online). In 
this case, the results of the J, H and K band fits were used to fix the values of the bulge parameters and the flux of 
the nucleus was estimated as Fnuc = Fgai — Ff,ui ■ 

MGG -2-40-4 — The bar and the disk are indistinguishable for r > 2" so the fitted component in the L-band is a 
mixture of both structures. The bulge parameters in this band were constrained to values found in the NIR to reliably 
separate it from the mixed bar-disk component. The M band profile was fitted with the bulge effective radius and 
Sersic's index constrained to the values found in the NIR. 

IRAS I548O-O344 — Barely no restriction was needed to obtain a good fit in NIR bands. In the L band only a nucleus 
and a bulge were fitted, and the index of the Sersic's law was limited by the values found in the NIR bands. 

IRAS 19254-7245 — This is an interacting system whose two nuclei are too close to fit reliably. Hence, no SED was 
obtained. 

NGG6810 — To obtain a good fit for this galaxy bumps in the radial profile produced by two large star forming regions 
were not considered. In the L, M and N bands, the seeing was better in the observations of the PSF star than in the 
observations of the galaxy. It is impossible for the bulge to account for the difference at r < 0.5". Fixed values for the 
bulge parameters were used, similar to the values found at shorter wavelengths, and the nuclear flux was estimated as 
Fnuc = Fgai — Fbui ■ The extended emission of this galaxy in the N band is an exceptional case. 

NGC6890 — To fit this galaxy the bumps at r w 15" and r k, 30" were excluded from the fitting process because they 
correspond to the ending of the flocculent arms. In the L band the signal beyond r k, 4" corresponds to sky emission. 
The galaxy in the M and N bands is consistent with being point-like, so no fit is performed. 

IG5063 — This galaxy is classified as SAO, but the best fit uses only bulge and nucleus components (no disk), which 
get stronger from shorter to longer wavelengths. In the L band we restricted the values of the parameters of the bulge 
to be close to the values found in the J, H and K fits. In the M band a change in seeing is evident, causing the PSF 
star to be sharper and narrower than the galaxy. The bulge parameters were frozen and only the nucleus was allowed 
to vary. The galaxy is point-like in the N band, so no fit is needed. 
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MRK897 — In the L band the object fell too close to the edge of the image and it was not possible to obtain the 
profile necessary to accurately determine the bulge contribution, so the bulge parameters were restricted to have values 
similar to the J H and K ones. 

NGC7130 — The only restriction needed to fit the NIR bands was that the disk scale length was larger than the 
bulge effective radius. In the L band the seeing apparently varied between observations of the galaxy and the PSF 
star. In order to obtain a reliable nuclear flux, the bulge component was fitted between 2" and 4" and the parameters 
constrained to be as close as possible to the J, H and K values. However, the fit to the galaxy was not adequate 
(mainly in the innermost region) , so the nuclear fiux is estimated as Fnuc = Fgai — Fhui ■ In the M band, the galaxy 
looks almost point-like so no fit is needed. 

MCG -3-58-7 — The disk size in the K band image was constrained to be larger than the bulge size. In the L band, 
only a nucleus and a bulge were included, and the index of the bulge was constrained to have a value similar to those 
obtained at shorter wavelengths. 

NGC7496 — While fitting this galaxy, the K band bulge index had to be restricted by the values obtained in the J 
and H bands. Also, the disk scale length was restricted to be larger than the bulge effective radius. In the L and M 
bands only a bulge and a nucleus were included. In L the bulge size was restricted by the values found in the NIR 
bands; in M the size and index of the bulge component were constrained by values found at shorter wavelengths. 

NGC7582 — The NIR images show a prominent galaxy, orientated nearly face-on, and two very weak arms with a 
different apparent inchna tion angle. In the literature this system has been characterized by an inclination angle of ~ 60 
degrees ([Morris et al.lll985t) , corresponding to the disk trace d by the weak arms. The galaxy has a strong star- forming 
region near the nucleus, which is also observable in X-rays ()Bianchi et al.ll2007D . so the region near r w 2.5" was not 
included in the fitting process. 

NGC7590 — This galaxy was observed twice in the NIR, once in service mode and once in visitor mode, with a different 
detector each time. The two fits gave consistent results. The point source detected in the NIR bands contributes at 
most 0.3% to the total fiux. In the L and M bands, no point source was detected, again suggesting a very weak AGN. 
After subtracting the star forming component (see Paper II) a Spitzer-IRS spectrum of NGC 7590 showed only a weak 
AGN continuum. 

CGCG 381-051 — This object was not detected in the L and M bands. In the J, H and K bands the fit did not include 
the region between r « 15" and r k, 20", which correspond to the end of the galactic arms. The H band results were 
used to contains the bulge index in the K band and the bulge size in the J band. Although the bulge is small it can 
be distinguished from the bar in all images. 



